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We demonstrate that beyond the universal regime corre- 
lators of quantum spectral determinants A(e) = det(e — H) 
of chaotic systems, defined through an averaging over a wide 
energy interval, are determined by the underlying classical dy- 
namics through the spectral determinant 1/Z{z) = det(2 — £), 
where e~^^ is the Perron-Frobenius operator. Application of 
these results to the Riemann zeta function, allows us to con- 
jecture new relations satisfied by this function. 

One of the central themes in the study of the quantal 
behavior of systems which have chaotic classical coun- 
terparts, is the quest for a correspondence between their 
quantal and classical properties. Examples of such cor- 
respondences are the random matrix theory (RMT) |^] 
behavior of chaotic systems over small energy intervals 

, and the "scar" phenomenon in which wave functions 
of chaotic systems are peaked in the vicinity of short iso- 
lated unstable classical orbits |^ . The universal RMT be- 
havior of chaotic systems is associated with the long time 
dynamics of the chaotic system, when classical probabil- 
ity distributions have already decayed into the ergodic 
stationary state. In the field theoretic approach to disor- 
dered metals Q RMT emerges as the zero mode contri- 
bution associated with the ergodic state. Scars, on the 
other hand, manifest individual non-universal properties 
of chaotic systems associated with the short time behav- 
ior of the classical system. 

The correspondence between the quantum spectral 
statistics and the classical behavior of chaotic systems 
is understood, usually, on a qualitative level. The pur- 
pose of this letter is to quantify this correspondence by 
formulating relations between spectral determinants as- 
sociated with both the quantal and the classical dynam- 
ics of chaotic systems. More precisely, we will show that 
correlators of spectral determinants of a given quantum 
system, defined through an averaging over a large energy 
interval, can be expressed in terms of spectral determi- 
nants of the classical counterpart. 

Correlators of spectral determinants of quantum 
chaotic systems have been previously discussed in the 
literature The results, however, were limited to the 
universal regime (zero mode approximation) where RMT 
applies. Therefore, they are incapable of providing any 
information regarding individual features of the under- 
lying classical dynamics. To this end it is crucial to 
go beyond the zero mode approximation, to encompass 
non-universal properties of the system. One possibility is 
to use periodic orbit theory. The non-universal proper- 
ties, in this approach, are related to short classical orbits. 



This approach has drawbacks which will be discussed be- 
low. 

In this paper we will employ an alternative semiclas- 
sical approach based on a recently developed supersym- 
metric non-linear a model for chaotic systems |6| which 
is free of these drawbacks. The basic classical ingredi- 
ents of this field theory are modes of the time evolution 
of the underlying classical system. These are the Perron- 
Frobenius modes, by which a disturbance in the classical 
probability density of chaotic system relaxes into the er- 
godic distribution [Q. They are analogous to the diffu- 
sion modes in the case of disordered systems Q . For time 
t > 0, the /i-th Perron-Frobenius mode decays as 6""*"^ * 
where 5R7^i > 0. 7^ are known as the Ruelle resonances. 
The dynamical zeta function 1/ Z[z) is the classical spec- 
tral determinant defined as 



l/Z{z)^JlA^{z-^^), 



(1) 



where are regularization factors introduced to make 
the product converge. This function, for systems which 
are ergodic on the energy shell, has a simple zero at the 
origin, 1/Z(0) = 0. This zero is associated with the 
ergodic state to which probability densities evolve, and 
manifests the conservation of probability. Thus Z{z), in 
the limit z — )■ 0, takes the form 



lim Z{z) = -— + V 

z^O TTZ 



(2) 



where U and V are real constants. 

The quantum spectral determinant which we consider 
here is associated with the spectrum of the Hamiltonian 
H. It is a real function of the energy e whose zeros coin- 
cide with the eigenvalues e„ of H, 



A{e)^l[AUe-er, 



(3) 



where are regularization factors. 

The correlators of the quantum spectral determinants 
will be defined with respect to the energy averaging: 
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where eo and N denote, respectively, the center and the 
width of the energy band over which averaging is per- 
formed. For convenience we shall work with dimension- 
less quantities, i.e. all energies will be measured in units 
of the mean level spacing, and h = 1. Let 71 be the 
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first non- vanishing Perron- Frobenius eigenvalue. We will 
assume that 3fi7i <C A'' <C eo- Under these conditions, 
results are independent of N, and the classical dynamics 
is essentially the same over the band. 

We consider here chaotic systems with no discrete sym- 
metries belonging to the unitary ensemble (broken time 
reversal symmetry). The formulae which will be derived 
and discussed in the following are: 



A(e + s)A(e-s) 
WW) 
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^{Z{2is)e-'^-"''} , (5) 
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A(e -I- s)A(e - s) 
A{e + uj+)A{e-uj+) 
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which apply if s, a; ^ N. When s, w ^ 3?7i, Z{z) can be 
approximated by (^) and the RMT results ||^ are recov- 
ered. 

An example different from RMT in which these rela- 
tions are exact is the one dimensional harmonic oscillator 
H ~ {p^ +q^) /2. This system is not chaotic, nevertheless, 
it is ergodic and its spectrum e„ = n -I- 1/2 is rigid. The 
quantum spectral determinant at high energies is simply 
A(e) = cos(7re), while the classical dynamical zeta func- 
tion is l/z{z) = 1 - e^'"^ with U = V ^1/2. Expanding 
l/cos(e+) = 2 ^^g(— l)"e*^*(2"+-^' and averaging over 
the energy it is straightforward to check all three rela- 
tions above are satisfied exactly. Notice that in contrast 
with chaotic systems, in this system all zeros of 1/Z{z) 
are purely imaginary. 

The derivation of the spectral determinant relations for 
general chaotic systems consists of two elements: One is 
the field theoretic framework of a semiclassical approxi- 
mation in which the classical ingredients are the Perron- 
Frobenius modes Q , and the second is a method of eval- 
uating non-perturbative parts of correlators 

In order to sketch the derivation, we consider relation 
(0) first. Our starting point is a formula for the ratio of 
four spectral determinants. Introducing the eight compo- 
nent field = (V'fl>fi,Xi?,Xfl>A,V'A>XA,XA), where 
Tp and X denote commuting and anti-commuting variables 
respectively, while RIA designate retarded/ advanced 
components, we have 



A(£ + s+)A(e-.s+) 
A(e + cj+)A(e-tj+) ' 

I'^'expi / dq^\(i)K 



(8) 
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where q is spatial coordinate vector, 
A = diag(l, 1,1, 1,-1, -1,-1,-1), and F = {F,-F) 



with F = diag(a;"'", w"*", s+j 5+). From here we proceed 
along the route discussed in Ref. Namely, an energy 
averaging (^ is performed which induces a quartic inter- 
action among the vp fields. Then, this interaction is de- 
coupled by means of the Hubbard-Stratonovich transfor- 
mation, and the integration over vp is performed. Finally, 
a steepest descent approximation, in which the width of 
the band N plays the role of the large parameter, is used 
to integrate out the massive modes. The result is of the 
form of an integral over the remaining Goldston (mass- 
less) modes. 



)A(e-s^ 
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(9) 
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where Q = T^^AT, and T = T(x||) belongs to the coset 
space C/(l,l/2)/[J7(l/l)(g);7(l/l)] (lO|. Here X|| denotes 
a vector of phase space variables which lies on the energy 
shell eo = -ff(x), and normalized such that J dx\\ = 1. C 
is the generator of time evolution of the classical dynam- 
ics defined by the Poisson brackets. 



C- = {H, ■}. 



(11) 



The eigenvalues of C arising from the regularization of the 
functional integral (^ constitute the Perron-Frobenius 
spectrum {7^}. 

The large frequency assymptotics of the above in- 
tegral can be calculated by stationary phase following 
Ref. 1^. Perturbation theory corresponds to the inte- 
gration of small fluctuations of T around unity, so that 
T = (1 — iP)^^ where P is small, and 



A{1 + iP){l - iPy 
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However, Q = A is not the only stationary point on the 
coset space C/(l, l/2)/[f7(l/l)(8);7(l/l)]. For the unitary 
ensemble there exists a second stationary point which is 
g = fcA where k = diag(l, 1, -1, -1, 1, 1, -1, -1) §. In 
our case it amounts merely to a change of sign in the 
Fermion block, namely to the substitution of s — > —s 
in the final results obtained from the first stationary 
point. It is therefore sufficient to consider the fluctua- 
tions around Q = A. Expanding the action to the leading 
quadratic order in B and B we have 

5*6// ~ i27r(tj+ - s) -I- (13) 

+ J dxi\Str FBB + FBB + iB£B . 

Expressing -B(x|| ) and i?(x|| ) in terms of the eigenbasis of 
the left (p^ and right (f>^ Perron-Frobenius eigenfunctions, 
as B(x||) = Ep-S^'/'m(x||) and i?(x||) = J2„ B^4>^{k\\), 
we obtain 
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-(7m -is- it^+)Str(B^B^)} + i2tt{uj+ - s). (14) 

We remark, here, that the Perron-Frobenius eigen- 
modes he outside the Hilbert space where B{yi\\) are 
defined. Nevertheless, they form a biorthonormal set, 
/ (ix||0y(x||)(/)^(x||) = 6v,^i^ which spans this space 

Integrating over B^, and taking into account also the 
contribution of the second stationary point Q = fcA we 
obtain (|^). Our perturbation theory can be justified only 
in the large frequency assymptotics s,uu 3> 1, neverthe- 
less, comparing it with the exact universal expression at 
the low frequency limit |^ and with the result of renor- 
malization group treatment |ll| , we conclude that (Q) is 
correct for arbitrary values of s and uo. 

Relations (||) and (^) can be similarly obtained by 
performing the energy averaging only in the Fermion- 
Fcrmion or in the Boson-Boson blocks respectively. The 
normalization constants of the resulting expressions are 
determined by considering the limits s,u}^ — > 0. 

It is instructive to show how one might use periodic or- 
bit theory in order to argue that the spectral determinant 
relations are satisfied. One can express both spectral 
determinants, A(e) (in the semiclassical approximation) 
and 1/Z{z), in terms of the classical periodic orbits of the 
system. For concreteness we consider two dimensional 
systems. The spectral determinant (^), in the semiclas- 
sical approximation, can be written as W3 



A(e)-e— ^(^)Cs(e), 



(15) 



were N{e) is the mean level staircase, and Cs(e) is the 
Selberg zeta function. 



p k=0 \ ' P 
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Here p labels the primitive periodic orbits of the system 
characterized by action Sp{e), Maslov phase i^p, and the 
eigenvalue of the monodromy matrix Ap with absolute 
value larger than one. 

The dynamical zeta function (|l|) has an exact repre- 
sentation in terms of a product over periodic orbits 



p fc=0 



|Ap|A^ 



fc+i 



(17) 



where Tp — dSp{e)/de is the period of the p-th orbit. 

Both periodic orbit products, (|l^) and (p7|), suffer from 
convergence problems in the regions of interest, neverthe- 
less they imply a particular form of regularization since 
they can be understood as analytically continued from 
domains of the complex plain where the products con- 
verge. Algebraic manipulation with their forms should 
be therefore understood as either performed using re- 
summed formulae, or analytically continued from regions 



of complex energies (e and z) where the products con- 
verge. 

To calculate (A(e + s)A(e — s)) one might try to ex- 
press (|l^) as a Dirichlet sum over pseudo-orbits , and 
then keep the diagonal part of the resulting double sum 
in A(e + s)A(e — s). However such calculation yields a 
wrong result. The reason is that the pseudo-orbits in 
the Dirichlet representation of the the Selberg zeta func- 
tion are not independent. They are related by the fact 
that A(e) for real e is a real function while its divergent 
periodic orbit representation is not manifestly real. In 
fact, a resummation of the divergent tail the Dirichlet 
sum reproduces the complex conjugate of the head, and 
the resulting expression for the spectral determinant has 
a Riemann-Siegel form in which only a finite number 
of orbits contribute effectively Using such a for- 

mula for A(e) one can show that to the leading order 
in l/|Ap| (||) is reproduced. In the case of inverse spec- 
tral determinants one should interpret l/Cs{E + lo^) as 
the usual periodic orbit product ( p^ ) evaluated at lo with 
sufficiently large imaginary part. Applying the diagonal 
approximation with the above interpretation of the spec- 
tral determinants at the numerator and the denominator 
and using |Aj,| ^ 1, one can obtain ^ and (^. 

It is worth pointing out the differences between 
the field theoretic and the traditional semiclassical ap- 
proaches. The starting point of the periodic orbit theory 
is a semiclassical approximation for the quantum spec- 
tral determinants given in terms of periodic orbits (p^). 
Energy averaging is then performed using the diagonal 
approximation. The problem is that the diagonal approx- 
imation is uncontrolled. Moreover it leads to wrong re- 
sults whenever the diverging periodic orbit sums involved 
in the calculation are interpreted incorrectly. The results 
of this approach also appear to be only a leading order 
in the instability of the orbits namely when |Ap| ^ 1. In 
the field theoretic approach the starting point is an ex- 
act quantum mechanical expression for the spectral de- 
terminants. The energy averaging is performed exactly, 
and only then the semiclassical approximation is applied. 
Here the approximations at each stage of the calculation 
are controlled and large instability, |Ap| ^ 1, is not re- 
quired. 

Our results can be used in order to reveal new mathe- 
matical relations satisfied by the Riemann zeta function. 



Cix) 



(18) 



where the product is over all prime numbers p. In view of 
the Riemann hypothesis, it is believed that there exists 
a chaotic Hamiltonian whose eigenenergies coincide with 
the non-trivial zeros of C{x) on the critical line — 
1/2 0. Thus C(l/2 - ie) plays the role of the Selberg 
zeta function for this system. The spectral determinant 
(which is a real function for real energies) is given by |lj] 



(19) 
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where N{e) = ^ (in ^ - l) + | + O (7) is the mean 
counting function, which gives the mean number of zeros 
below e. The mean density of zeros d{e) = dN{e)/d€ ^ 
ln(e)/27r, thus, increases logarithmically with the energy. 
Here, for convenience, we choose not to normalize the 
energy by the mean level spacing. 

As we have shown, correlators of quantum spectral de- 
terminants are related to the classical spectral determi- 
nant. To apply our results it is therefore necessary to 
determine the function 1/Zr{z) which plays the role of 
the dynamical zeta function for the "Riemann system" . 
This can be deduced from the exact periodic orbit for- 
mula ITll, 



satisfied by C,{s). Our relations (||) (||) and ^ can be 
also straightforwardly generalized to higher order corre- 
lators. 
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l/Z{z) ^ exp. 



EE-|det(M;-/)| 
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where Mp is the monodromy matrix of the p-th orbit 
repeated r times [for two dimensional systems Eq. ( |20| 
reduces to (|l^)]. From, Tp = logp and | det(M^ - /)| = 
grinp rra follows that 



Zfl(z) = C(l-z), 



(21) 



Thus the inverse of the dynamical zeta function is itself 
the Riemann zeta function. It has a single simple zero at 
the origin, and since lim^^o C(l — z) = —1/ z + C , U = tt 
and V = C where C — 0.577215 • • • is Euler's constant. 
Assuming the existence of an underlying chaotic "Rie- 
mann system" , and using ( p^ ) and (|2^) one can apply 
®; ® ^-nd (0)i to the Riemann zeta function . 

It is worth mentioning that the periodic orbit ap- 
proach, namely using the Dirichlet representation, ^(x) — 
^^j^n"^, the fact that Ai^(e) is real, and diag- 
onal approximation, yields (Afl(e -I- s)Afl(e — s)) — 
3fi|Zfl(2is)e-*2'rrf(eo)s| ^i^icij ig analogous to (|), and 

exact to any order in l/|Ap|. In this respect the diagonal 
approximation for the Riemann zeta function is compat- 
ible with the field theoretic calculation. 

Relation (^ for the Riemann zeta function can be writ- 
ten in a particularly simple form. Substituting ( |T9| ) and 
(|l|) in @ leads to 



= C(l - i2LU+). 

(22) 
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z(e 



To summarize, we have shown that correlators of quan- 
tum spectral determinants of chaotic systems are ex- 
pressed in terms of the classical spectral determinant as- 
sociated with the Perron-Frobenius spectrum. On en- 
ergy scales smaller than 3?7i where 71 is the first non- 
vanishing Perron-Frobenius eigenvalue, the universal re- 
sults are reproduced, while over larger scale the behavior 
is determined by properties of the dynamical zeta func- 
tion. The results, applied to the Riemann zeta function, 
enabled us to conjecture the existence of new relations 
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